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SUMMARY 


A s t iid v ol the available methods for transonic airfoil and wing design 
indicates t lt.it the most powerful technique Is the numerical optimization pro- 
cedure. However, the computer time for tills method Is relatively large 
because of the amount of computation required in the searches during optimiza- 
tion. The optimization method requires that "base" and "calibration" solu- 
tions be computed to determine a "minimum drag" direction. The design space 
Is then computationally searched in this direction; It is these searches that 
dominate the computation time, A recent similarity theory allows certain 
transonic flows to be calculated rapidly from the base and calibration solu- 
tions. In this report the application of the similarity theory to design prob- 
lems is examined with the object of at least partially eliminating the costly 
searches of the design optimization method. An example of an airfoil design 
Is presented. 


INTRODUCTION 


New interest in reducing the fuel consumption of existing and future alr- 
craft has been prompted by the world-wide fuel shortage. A logical place to 
start the search for lower fuel consumption Is to reduce drag, especially at 
the transonic speeds at which most aircraft cruise. The most common descrip- 
tion ot a transonic flow is when there is a supersonic "bubble" totally 
embedded in a subsonic flow. The superronic bubble may be terminated by a 
shock wave producing wave drag or, in certain circumstances, may return to 
subsonic conditions through an Isentroplc compression with no wave drag. This 
wave drag is associated with the entropy change across the shock. Soon after 
a shock wave appears in the flow, the drag will increase rapidly with increas- 
ing tree-stream Mach number, leading to the definition of a "drag rise Mach 
number," which is defined as the free-stream Mach number at which this rapid 
drag rise begins. One of the main objects of designing a wing for transonic 
speeds is to obtain as high a "drag rise Mach number" as possible, subject 
to certain constraints. The obvious way to reduce the wave drag, at least in 
a two-dimensional study, is to use a supercritical shock-free airfoil section 
where there is no shock wave and, consequently, no wave drag. However, these 
shock-free airfoils may have undesirable off-design characteristics, such as 
strong shock waves suddenly appearing when the Mach number is perturbed 
slightly from its design value. Thus, an important constraint in the design 
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of alrtollH 1h that there Hhuuld be good off-design behavior. In wlnR design 
(i.e., three-dimensional flows) the deHlRti can be altered by span wise changes 
(e.g., sweep or twist); and this can complicate the design procedure. 

Methods for calculating transonic flow characteristics around given air- 
foils were first derived bv Murman and Cole (ref. 1). with a later correction 
by Murman (ref. 2), using the concept of type-dependent finite differences 
applied to the transonic small disturbance equation. Later work resulted in 
methods using the full potential equation such as that by Jameson (ref. 3). 
Three-dimensional extensions of the small disturbance methods were made by 
liallhaus, Bailey, and Frick (ref. 4), and extensions of the full potential 
equation methods were made by Jameson and Caughey (ref. 5). All these meth- 
ods are for lsentropic flow, which, strictly speaking, cannot produce a wave 
drag that (as noted above) would be a consequence of an entropy change. How- 
ever, lsentropic flow does produce a momentum deficit across the shock wave 
that can le thought of as a wave drag (ref. 6), an assumption that seems to 
be corroborated by experiment. These calculation methods must be differenced 
(ref. 2) In conservative form otherwise the shock location and strength muy 
be incorrect. 

The earliest attempt to design transonic airfoil sections used the hodo- 
graph transform.it ion, in which the velocity components are the Independent 
variables rather than the usual geometric variables (x,y). The hodograph 
transformation leads to a linear equation for transonic I low rather than the 
usual nonlinear equation in the physical coordinate system. Best known of 
these methods are those of Nlcuwland (ref. 7) and Bauer, Cambodian, and Korn 
(ref. 8). These methods require complex mappings and transformations and can 
be difficult to use. In addition, neither the design pressure distribution 
or the design Mach number Is known in advance. Only shock-free designs are 
possible, thus precluding the Inclusion of off-design criteria into the 
design process. Also, It is difficult to apply constraints, such as maintain- 
ing a specified minimum lift coefficient. 

A second type of design method is the reversal o! the direct finite dif- 
ference procedures, with the pressure distribution specified at a given Mach 
number and the corresponding airfoil shape obtained. Examples of this type 
of procedure are given by Steger and Kllneberg (ref. 9) ar.d by Carlson 
(ref. 10). The main problems associated with this line of attack are: 

1. The specified pressure may not produce a realistic airfoil (e.g., 
nonclosure of the trailing edge). 

2. Constraints on the lift coefficient or section thickness are diffi- 
cult to implement. 

J. Off-design criteria cannot be incorporated into the design process. 

A third method of transonic wing design is the numerical optimization 
procedure developed by Hicks (refs. 11 and 12) and his coworkers. The numer- 
ical optimization method s*-eks to minimize some specific parameter (e.R., the 
drag coefficient, for a set ot design variables describing the airfoil 





geometry, while satisfying .1 number of specified constraint*. The constraints 
may b** aerodynamic (e.g., 1 If t -coef f Ic lent or off-design criteria) or geo- 
metric le.g., airfoil thickness or volume). The design method uaea an aero- 
dynamic analysis program coupled with a numerical optimization program; the 
design process is briefly outlined below. 

A base airfoil profile is chosen and is perturbed by the use ot shape 
functions which control the final profile. The coefficients of theae shape 
functions are the design variables. The optimization scheme perturbs each of 
these coefficients in turn, returning to the aerodynamic analysis program for 
the evaluation of the drag coefficient after each perturbation. After the 
change ot Cp lor each change oi design variable lias been noted, the optimi- 
zation program computes the gradient of Cp (VCn) with respect to each design 
variable. The optimization program then increments the coefficients one to 
four times in the -VCy direction, searching for a minimum value of that 

satisties the constraints. At each increment, the aerodynamic analysis pro- 
gram is used to compute the drag and aerodynamic characteristics and any off- 
design constraints. If a constraint is reached or the drag increases (due to 
nonlinearity), then a new gradient is computed and the process repeated. The 
computing time of this optimization procedure Is fairly large. The bulk of 
the computing time, however. Is used In searching for the minimum drag, when 
the analysis program has to be used at each incremental step. The inclusion 
of off-design criteria in the constraints require an additional calculation 
of the aerodynamic characteristics at the specified off-design Mach number. 

An outcome of this design effort (refs. 11 and 12) is the determination of 
useful airfoil shape functions that, when added to the base airfoil, permit 
a large class of airfoil contours; these shape functions always give closure 
of the airfoil. Extensions ot the numerical optimization procedure to finite 
wings has been reported bv Hicks and Hcnne (ref. 13). In this three- 
dimensional procedure, five shape functions are used at two span-wise sta- 
tions, together with an angle of attack variation, giving a total of eleven 
design parameters. While angle ot attack and wing twist are computationally 
efficient to use in the design code, other design criteria (such as sweep, 
aspect ratio, and taper) require a great expenditure ot computing power. This 
is because these design variables affect the location of the wing relative to 
the finite difference mesh (which is usually sheared so that the grid lines 
coincide with the leading and trailing edges), and modifications to the mesh 
are desirable for each change. Another application of this design procedure 
is by Haney, Waggoner, and Ballhaus (ref. 14). 

It would seem from the above discussion that the most powerful of the 
available design methods is the numerical optimization technique, since aero- 
dynamic constraints such as lift coefficient, off-design criteria, and geo- 
metric constraints (such as wing thickness) can be Incorporated into the 
design process. However, the computing requirements can be considerable, 
especially for three-dimensional designs. The bulk of the computing time is 
used in the searches where the aerodynamic analysis program must be used at 
each step. If some way of reducing the time for the searches could be found, 
then the computer requirements would be less formidable. 
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Recently, a similarity theory has been derived bv Nixon (refs. 15 and 16) 
by which a range of transonic flow solutions can be found by simple algebra, 
provided two solutions, a base solution and a calibration solution, are known. 
The analysis is based on a perturbation expansion of the transonic equations, 
which leads to a linear oerturbatlon equation. The nonlinear phenomena of 
shock movement due to a perturbation Is treated bv using a strained coordinate 
system in which the shock location is invariant. Essentially, this device 
treats the nonlinear flow changes due to a perturbation in some flow charac- 
teristic, such as Mach number, by a nonlinear combination of two linear prob- 
lems that can be solved In sequence. Since these equations are linear, the 
principle of superposition applies; and hence, the effect of several pertur- 
bations can be considered at once. The nonlinearity appears only in the last 
step of the procedure, which is the nonlinear combination of the linear prob- 
lems. Generally, if N parameters are perturbed, then the procedure requires 
N + 1 solutions composed of one base solution and N calibration solutions. 
Once the*e solutions are known, any related "nearby" solution can be obtained 
rapidly. Using a CDC-7600 computer for a three-dimensional, two-parameter 
example, eight cases can be computed in 0.38 s of CPU time. The theory lias 
been applied to three-dimensional wings (ref. lb) with multiple shock waves. 
The main restriction of the theory is that shock waves cannot be generated or 
destroyed in the perturbation; although in principle a shock of zero strength 
can be treated, that is, a shock-free supercritical airfoil. 

Since the computation time required to calculate the "nearby" flow char- 
acteristics is so small, it seems that this similarity method would be 
extremely advantageous in the numerical optimization design procedure. The 
costly searches previously carried out by reversion to t lie aerodynamic analy- 
sis program can now be done using the inexpensive similarity theory. Also, 
since the similarity theory gives the "nearby" solutions as analytic rather 
than numerical functions of the parameters, it is possible that the optimiza- 
tion procedure itself could be improved by a study of the analytic dependence 
of the design criteria on these parameters. Furthermore, off-design con- 
straints in the transonic regime, e.g., drag at a slightly lower Mach number, 
can be easily computed since a change In the Mach number is just another per- 
turbation parameter. This avoids the need to calculate off-design character- 
istics at each point in the design loop using the aerodynamic program. 

Another example is the effect of wing sweep angle in three-dimensional flows; 
which, again, is Just another perturbation parameter, the effect of which need 
only be computed once. 

This report begins to investigate the applications of the similarity 
theory to design problems. Only two-dimensional flows are considered, and a 
total of five parameters are used to characterize the shape functions. At 
this stage an optimization procedure is not used, since the basic aim is to 
establish the validity of the design applications of the theory. The main 
objective in the present design study is to reduce the wave drag coefficient 
of an airfoil and to deduce the necessary ground rules for this design objec- 
tive. Applications of these rules to the design of an airfoil section are 
presented, and the design pressure distribution agrees satisfactorily with a 
direct calculation. 
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HAS 1C PRINCIPLES 


It an nirfoll section in specified bv several shape function* with 
unknown amp lit tides , the problem under consideration is to choose these .impll- 
tudes such that Home design criteria 1h met. General iv, in the course 
choosing these design parameters, anv shock waves in the flow will change 
location, thus invalidating the usual type of perturbation theory. This is 
because tho pressure changes in the region traversed by the shock wave are not 
small, even for a small change in shock locations. A means of treating this 
kind of perturbation is given by Nixon (ret. 15), using a strained coordinate 
system in which the shock location is invariant with changes in the perturba- 
tion variable. Using this technique, the nonlinear transonic problem is split 
into two linear problems which can be solved in sequence. Because ol this 
linearity, complex solutions can be constructed from simpler solutions bv 
superposition. In a design problem the effect of each design parameter can 
therefore be calculated separately and the principle of superposition used to 
construct a solution giving the effect of any combination of the design 
parameters. 

The analysis given by Nixon (ref. 1*>> concerns only the transonU small 
disturbance equation and since design methods are most at urate if the full 
potential equation is used, the strained coordinate technique must be devel- 
oped for such an equation. 

It (x,y) is a Cartesian coordinate system, nondlmenslonal ized with 
respect to the airfoil chord and with x aligned with the airfoil chord, then 
the full potential equation is 

(a’ - u ; )$ - 2uv$ + (a*’ - v 3 )i • 0 (1) 

xx xy yy 

where u and v are the velocity components in the x and y directions, 
respectively, nondlmenslonal lzed with respect to tin tree-stream velocity 
‘L. • are given by 


u 


. It 

3x 


v 


3 £ 

“ 3y 


( 2 ) 


where $ is the velocity potential; a Is the speed ot sound, again nondl- 
mensionalized with respect to q,„ and given by 


.i 


2 



- v ) 


(3) 


M», is the tree-stream Mach number. The boundary condition on the airfoil 
surface is given by 
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v(x,y 8 ) 


( 4 ) 


u(x,y M ) 


y 8 (x) 
X 


where y ■ v lx) defines the surface of the airfoil. An infinite distance 
upstream of the airfoil 


u ■ q^coxa 

v • q xlna (5) 


where a is the angle of attJ*' ! . 

The pressure velocity relation is given by 


Y /y_j 

C p (x,y) ■ ^ Mb, 7 £ L ♦ (1 - u 2 - v 2 )] (6) 

It is now proposed to change the airfoil boundary by some small amount char- 
acterized by » . Thus, the new airfoil ix defined by 

y • ? H (x) * y B (x) + C Y ^(*) (7) 

aid the corresponding boundary condition lx 


v(x,y s ) 

“(x.v ) 


V x) 


+ 



(X) 


(8) 


The question is, given the solution of the problem defined by equa- 
tions (1) and (4), find the solution ol the problem defined by equations (1) 
and (8). As in reference lb, it is assumed that any shock waves in the flow 
are normal to the x-axis and that there is only one shock wave on each airfoil 
surface. The strained coordinate x' is tb n Introduced where 


x* + e 6 x x.(x') 

S 1 


(9) 


where e6x^ is the shock movement and xi(x') is given in reference 15 by 


, n x' (1 - x' ) 
»i ( x ) ■ x ,{, : 

S 8 

XI (x' ) - 0 


0 c x' < 1 

x’ < 0 
x' > 1 


( 10 ) 
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wliero x; 1m the original stock location and tfix^ 1 m the shock movement. 
The variables f, u, etc., are then expanded as the Merles 


♦(x,y) - ♦ 0 <*'»y‘j+ (*' »y) ♦ . . . 

and 

u(x,y) - u o (x\y)£l - c«x r xi ^ (x' )1 ♦ cui(x'.y) . . . 
v(x,y) • v^x'.y) ♦ rvjtx'.y) . . . 


(ID 


( 12 ) 


a(x,y) - a o (x',y) ♦ cai(x'.y) . . . ( 13 ) 

The expression for a(x,y) can be found In terms of u, v, and fix , using 
equations (3) and (9). 8 

Substituting equations (9), (11), (12), and (13) into equations (1) and 
(8) and equating coefficients of c gives, to first order in c, the follow- 
ing two equations with associated boundary conditions 



u ‘ )♦ 
o o,, 
x'x' 


2u v 4> 
o o o , 

X v 



V 2 
o 



0 


V o (x, »>’s ) 

u o (x\y g ) 



(x 1 ) 


(14) 


(15) 


and 


( a * “ u • 2u v $ , + (a • - v * )♦, + 2(a a, - u u, )♦ 

o o lx x o o T lx'y o o lyy o 1 « 1 

“ 2(u v ♦ v u )♦ + 2 (a a, - v v, )♦ 

o 1 o 1 o mIi> o 1 o 1 oyy 


xx 


x y 


■ “2$ u Sx x, , - 2$ fix x, , (a - u 2 ) 
O x , x , O 8 x' s 1 x' O o 


X X 


- (a 2 - u 2 ) fi x x $ + 2u V fix X, * 

o O 8 1 , t O 1 0081,0, 

x'x x l x' x'y 


(16) 


v l (x ’* y 8 ) “ (x')u o (x',y 8 ) 9 - ( U1 - fiX B U o X 1 )y (x‘) 

X x’' x' 

+ , ,<„•) (O) 


X X 
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Provided the perturbation is not i function of M„, then equations (lb) 
and (17) do not contain the parameter c. Hence, if the variables , 6x H , 
etc., ire kno'm for one value of c, then the value of the variables for any 
other value of t can be found by simple proportion. Chus, 


* Uj(x' ,y) 
c Ax 

s 



( 18 ) 


The omission of the Much number variation from the analysis Is because the 
term ai contains M*,; hence equation (16) is not independent of c. How- 
ever, Mach number variations can be treated if it Is assumed that for such 
changes the transonic small disturbance equation Is a valid approximation. In 
this caue, the analysis of reference 16 can be used. 

Wlille the linear perturbation equation, equation (16), can be solved for 
U} , 6x h , etc., it is much more convenient to use the same te hnlque to solve 
both tiie base and calibration solutions. Equation (16) multiplied by c c 
represents the first order of magnitude in c the difference of two solu- 
tions to equation (1). Hence, expressions for ii) and Ax s in equation (13) 
can be found by ,1 suitable combination of two known nonlinear results. 

If, by some method, the solution to the base problem defined by equa- 
tions (14) and (15) is known and If the solution to some perturbed problem, 
characterized by some parameter * , is also known (the calibration solution), 

then the terms (c ui ) and (t Ax ) can be found as follows: 
o * os 

1. The change In the shock location c 0 Ax„ between the base and cali- 
bration solutions is easily tound by Inspection. 

2. If u^)(x,y) is the soli. t ion of the perturbed problem and if 
u (0)(x*,y) i 8 the solution of the base problem, then 


v, (x ' * y) 


u< 1 )(x,y) 



c Ax x. (x')l 
° S l x , J 


(19) 


where 


x ■ x' + e Ax x. (x* ) 
o s 1 


( 20 ) 


and xi(x') is given by equation (10). 

Having obtained Ax and Uj(x',y) ore can then obtain the final 
solutions 


u(x,y) 


u( 0 )(x’ 



cAx x, 
s 1 



♦ cuj (x' ,y) 


( 21 ) 
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and 


x • x' 4- c«x m x ) (x') (22) 

The other velocity component la h I mply given by 

vi(x,y) - 7~(vt l Ux.y) - v< 0 )(x,y)) (23) 

c o 

where v(°)(x,y) and v^’^(x,y) are the aolutlona of the bnae and calibration 
problems, respectively. The total velocity in the y direction 1 h then 
given by 


v(x,y) • v(°)(x,y) ♦ cvj (x,y) (24) 

The pressure coefficient can then be found from equations (6), (21), 

(22), and (24). 

Since the basic equation derived in the preceding sections are linear, 
the effect of more than one purameter can be obtained by superposition. 

Thus, for N parameters, equations (21), (22), and (24) can be generalized 
to give 

f N 1 N 

u(x,y) - u(°)(x',y) 1 - £ e i 6x u x , '» x ') ♦ Za c i u < (*'•*) (25) 

L l-i i x* J f-i 1 l i 

where 

N 

x • x' + 2 c . fix x j(x') (26) 

1-1 1 8 i 


v(x,y) - v (o) (x,y) ♦ £ c v. (x,y) 

l-l 1 i 


(27) 


The N parameters are denoted by c^(i"l,N), the change in shock locution 
due to the ith parameter change is c.6x 

1 8 i 

DESIGN APPLICATIONS 


In the direct calculations, the parameters e are specified and u(x,y) 
and v(x,y) on the airfoil surface are obtained. In a design application, 
the velocity u(x,y) on the airfoil surface can be specified at N - 1 
stations and the c found; the Nth equation for the c. is found by 

° N 1 

specifying the total shock movement £ c.6x , thus giving the relationship 

i-1 1 8 i 
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between the coordinate sy«.icms x and x'. Once u(x,y) on the .ilrfoll 
mir f .ue and the i are known, the v(x,y) can be obtained from the tangency 
boundary conditions. An alternative to apecifylnft u(x,y) la to uae eone form 
of optimization procedure. The parametera t ( may be changea in angle of 
attack, geometric ahape functlona, or, If thesmail disturbance equation la 
lined , changea In Mach number. 

Am shown in reference 16, the adequate calculation of the drag coeffi- 
cient Involves a flow field calculation around the ahock waves and is n cubic 
equation in the parameter c. This would lead to a complicated analysis In 
the design application. Therefore, for optimization purposes. It Is proposed 
that the possibility of minimizing the function | Cp - * - - t *| rather than Cp 
is Investigated, where Cp+ is the surface pressure Just ahead of the shock 
wave and C * Is the critical pressure. The applicability of this assumption 
rests on c{] being a monotonic function of |Cp+ - C p *| . In figure 1, a plot 
of Cp against |Cp + - C.,*t for 14 different direct calculations is shown, 
and It may be seen that within the limlts^of numerical accuracy Cp is 
indeed a monotonlc function of |Cp* - Cp | . 

An obvious design objective is to reduce the drag coefficient by reducing 
the shock strength to zero or nearly zeto. In theory, this should give some- 
thing close to a shock-free airfoil. In order to establish some ground rules 
for such a design objective, a test case using the transonic small disturbance 
equation with linear boundary conditions is considered. Tae base airfoil is 
a 102 biconvex section at zero angle of atlack and the object is to find the 
ftce-stream Mach number at which the shock strength Is nearly zero. The base 
Mach number is 0.828 and the calibration Mach number is 0.838. A shock fit- 
ting small disturbance code was used to compute these solutions. The design 
Mach number is 0.7905 and the result obtained by the similarity solutions is 
compared to a direct calculation in figure 2. The difference In shock loca- 
tions and pressure distributions ahead of the shock are probably due to the 
magnitude of the perturbation In Mach number being too large. However, the 
main discrepancy between the similarity and direct results is in the large 
supersonic expansion In the former behind the shock wave. This expansion is 
due to the postshock expansion behavior exhibited by the base and calibration 
solutions, which apparently does not scale when a shock-free limit is 
approached. This suggests that in order to adequately compute shock-free or 
nearly shock-free solutions, both buse and calibration solutions should fairly 
closely model the essential flow features of tie final design. For example, 
the rapid postshock expansion should not be large in the tase and calibration 
solutions. In figure 3, the pressure distribution close to the shock at 
M^, ■ 0.798 for the Bame biconvex airfoil is shown, but with a base Mach num- 
ber of 0.818 and a calibration Mach number of 0.808. It can be seen that the 
rapid postshock expansion has been eliminated. These results lead to the 
following design criteria for shock-free or nearlv shock-free designs. 

1. The base and calibration solutions should represent all of the essen- 
tial features of the final design. 

2. The flow must not accelerate supersonically behind a shock wave 
(including a zero strength shock) and must be supersonic Just ahead of the 
shock (l.e., the shock must be compressive). 
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1, The final design should be realistic, that Is, there should be no 
crossover of the upper and lower surfaces of the airfoil. 

4. For a shock-free design, It Is desirable that any discontinuities In 
the pressure gradient at the "shock" be at a minimum In order to have a 
smooth recompression. 

A final criterion Is that the perturbat Ions should not be too large. 
Generally, this Is Indicated when the strained coordinate overshoots the air- 
foil chord, that 1s, gives values of x that lie outside the airfoil. Con- 
sequently, a fifth condition Is as follows: 

5. Fur values of the strained coordinate x' on the airfoil chord, the 
coordinate x must alwuys lie on the chord line. 

These then are the criteria used in the preliminary tests of design 
applications of the similarity theory. 


EXAMPLE 


In order to test the above ideas, an airfoil Is designed using the Ideas 
of Hicks (re' 11 and 12). The base airfoil Is a laminar flow design. 

Pertv Nations of the form 

y ■ a ^ (sin nxbj) 3 1 ■ 1,5 (28) 

are used to modify the airfoil geometry. The values of b^ and a^ used 
in computing the calibration solutions are shown In table 1. Both base and 
calibration solutions were computed using a full potential equation code. 

The object is to reduce the shock strength and hence, the wave drag by 
noosing n^d • 1,5). The magnitude of the perturbations is limited to 
1-1/2 times that used In the base and calibration solutions, although the sign 
can change. This is effectively a constraint on section thickness. The 
free-stream Mach number is 0.74 and the angl • of attack is zero. The optimi- 
zation Is simply done by a search of the relevant range of parameters with 
12 steps in each rar:;e. This scheme 1h not by any means the best, but is easy 
to program and Is used here only to validate the theory. The magnitude of the 
parameters found by this procedure Is given In table 1. The resulting pres- 
sure distribution is shown in figure 4 and compared to a direct calculation. 

It can be seen that the predicted and direct calculations agree fairly well. 
The total computing time Is 3.7 s on a “DC-7600 computer, provided the base 
and calibration solutions are known. 

In this paper, only one design iteration Is considered; that is, only one 
set of base and calibration solutions is used. This serves to test the ideas; 
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and. In any cahv , or.? such computation may be sufficient. If the bane and 
calibration Holutioua do need to 1« recalculated, then the whole procedure 
can easily be repeated. 


CONCLUDING KKMARKS 


The similarity theory developed by Nixon (refs. 15 and 16) for direct 
calculations Is extended to design optimization problems. Some ground rules 
for the design of shock-free or nearly shock-free airfoils are deduced. A 
simple example is computed. The advantages of the scheme of incorporating 
the similarity theory into the numerical optimization design procedure are as 
follows : 


1. In the hitherto expensive searches, the computing time can be reduced 
considerably using the similarity theory. 

2. Off-design criteria can be easily and inexpensively Incorporated Into 
the design process. 

3. In three-dimensional applications, the difficult design parameters, 
such as wing sweep or taper, can be much more easily taken into account since 
only one calibration solution is required rather than the multiple finite dif- 
ference calculations required in the existing numerical optimization scheme. 

4. The analytic nature of the similarity theory may lead to improve- 
ments ia the optimization scheme itself. 

The author thanks T. L. Holst and R. M. Hicks of NASA Ames for their help 
and A. R. Seebass and N. J. Yu of the University of Arizona for the shock- 
fitting, small disturbance results. 
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TABLE 1. DESIGN PARAMETERS 


B 

b i 

(calibration) 

a j (design) 

i 

0.2314 

-0.003 

0.0034 

B 

0.3654 

0.0005 

0.0008 

B 

0.6603 

0.002 

-0.0013 


0.8681 

0.001 

0.0001 

5 

2.4094 

-0.002 

-0.0022 
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Figure 2 
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Figure 3.- Pressure distribution at the shock wave 10' biconvex airfoil; 

M -0.7*18. 
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Figure 4.- Airfoil design, a ■ 0°, M x ■ 0.74 (upper surface only). 
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